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SUMMARY 

Failures of many pressure vessels have been traced to surface cracks. Accurate stress 
analyses of these surface-cracked components are needed for reliable prediction of their 
crack-growth rates and fracture strengths. Because of the complexities of such problems, all 

investigators have used engineering estimates or approximate analytical methods to obtain the 
stress-intensity factors. 

A few three-dimensional stress analyses of semi-elliptical surface cracks in pressurized 
cylinders have been reported recently. However, these investigators considered only an 
internal surface crack with a crack-depth-to-crack-l.ngth ratio of 1/3 and a wall-thickness- 
to-vessel-radius ratio of 0.1. 

The purpose of this paper is to present stress-intensity factors for a wide range of 
semi-elliptical surface cracks on the inside of pressurized cylinders. The ratio of crack 
depth to crack length ranged from 0.2 to 1; the ratio of crack depth to wall thickness ranged 
from 0.2 to 0.8; and the ratio of wall thickness to vessel radius was 0.1 or 0.25. The 
stress-intensity factors were calculated by a three-dimensional finite-element method. The 
finite-element models employed singularity elements along the crack front and linear-strain 
elements elsewhere. The models had about 6500 degrees of freedom. The stress-intensity 
factors were evaluated from a nodal-force method. In this method, the nodal forces normal to 

the crack plane and ahead of the crack front were used to obtain the stress-intensity 
factors. 

An empirical equation for the stress-intensity factors was fitted to the results of the 
present analysis as a function of crack depth, crack length, wall thickness, and vessel 
radius. The equation applies over a wide range of configuration parameters and was within 
about 5 percent of the present results. 

The present results were compared to other analyses of internal surface cracks in cylin- 
ders. The surface-crack configuration had a crack-depth-to-crack-length ratio of 1/3 and a 
wall-thickness-to-vessel-radius ratio of 0.1. Results from the literature using a boundary- 
integral equation method vtere in good agreement (±2 percent) and those from a finite-element 
method were in fair agreement (±8 percent) with the present results. 

The stress-intensity factors and equations presented herein should be useful in corre- 
lating fatigue-crack-growth rates and in calculating fracture toughness for the surface crack 
in a pressurized cylinder. 


1. Introduction 


Failures of many pressure vessels have been traced to surface cracks. Accurate stress 
analyses of these surface-cracked components are needed for reliable prediction of their 
crack-growth rates and fracture strengths. Because of the complexities of such problems, all 
investigators have used engineering estimates or approximate analytical methods to obtain the 
stress-intensity factors. 

Some engineering estimates for the stress-intensity factors for surface cracks in pres- 
surized cylinders have been made by Underwood [1] and Kobayashi [2]. Their estimates did not 
include the effects of wall thickness. Recently, however, Kobayashi, Emery, Polvanich, and 
Love [3] have estimated stress-intensity factors for internal surface cracks that did include 
the effects of wall thickness. 

A few three-dimensional stress analyses of serai-elliptical surface cracks in pressurized 
cylinders have been reported recently. Atluri and Kathiresan [4] and McGowan and Raymund f5l 
used three-dimensional finite-element methods, while Heliot, Labbens, and Pellissier-Tanon 
[6] used the boundary-integral equation method, to obtain stress- intensity factor variations 
along the crack front for a limited range of configuration parameters. References [ 5 5 and 
[6] considered only an internal surface crack with a crack-depth-to-crack-length atio of 1/3 
and a vall-thickness-to-vessel-radius ratio of 0.1. 

The purpose of this paper Is to present Mode I stress-intensity factors, calculated by a 
three-dimensional finite— element method [7—9] , for a wide range of semi-elliptical surface 
cracks in pressurized cylinders. The cracks were located on the inside of the cylinders. 

The ratio of crack depth to wall thickness ranged from 0.2 to 0.8; the ratio of crack depth 
crack length ranged from 0.2 to 1; and the ratio of wall thickness to vessel radius was 
0.1 or 0.25. The stress-intensity factors were calculated by using a nodal-force method 
l 7—9] . An equation for the stress— intensity factors was also developed for a wide range of 
configuration parameters. The stress-intensity factor variations along the crack front are 
presented and, where possible, compared with other solutions from the literature. 

2. Symbols 

a depth of surface crack 

b half-length of pressurized cylinder 

c half-length of surface crack 

f ratio of boundary-correction factors (cylinder to flat plate) 

F boundary-correction factor for an internal surface crack 

th 

Gj ,Hj boundary-correction factors for j stress distribution on crack 
Kj stress-intensity factor (Mode I) 

p internal pressure on the cylinder 

Q shape factor for an elliptical crack 

R,R q inner and outer radii of cylinder 

t cylinder wall thickness 

x,y,z local Cartesian coordinates centered at crack mouth 
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X.Y.Z global Cartesian coordinates 

® angular measurement on cylinder 

a j applied stress on crack surfaces 

<*■ parametric angle of elliptical crack 

Three-Dimensional Analysis 

A surface crack 1„ an iut.rually pressurised cyriuder is sho™ i„ flg „ c . a. 
cyl uder of wall thickness r, lucerual radius R, ,„d l e „gth 2b, contains a seui-elUptical 
;ur ace crack of lc»g t b 2c and depth a on the inner surface of the cyiinder. The slross- 
-ntensity factors for the surface-crack configurations were obtained by using a three- 
dimensional finite— element analysis* 

3. 1 Finite-Element Idealization 

TWo types of elements (isoparametric . * singular [7]) were used in combination to model 
the cylindrical vessels. Figure 2(a) shows a typical finite-element model for an internal 
surface crack. The model, which employed nearly 6500 degrees of freedom, idealized one- 
eighth of the vessel (0 < 0 < 90 degrees and 0 < X < b) . The isoparametric (linear-strain, 
eight-noded hexahedron) elements were used everywhere except near the crack front, where 
eight singularity elements, each in the shape of a pentahedron, were used [8]. The singu- 
larity elements produced a square-root singularity in stress and strain at the crack front. 

A typical finite-element pattern along the crack plane is shown in figure 2(b). The finite- 
element models in the neighborhood of the crack were identical to those used in refer- 
ences [8] and [9] for surface cracks in flat plates, except that the models were curved to 
conform to the desired cylindrical shape. The vessel half-length-to-radius ratio, b/R, 
ranged from 1 to 5, the b/c ratio ranged from 10 to 50, and the t/R ratio was 0.1 or 
0.25. Further details on modeling and formulation of the two types of elements used are 
given in references [7-9] and are not repeated here. 

3- 2 Boundary Conditions and Applied Loading 

Symmetry boundary conditions were applied on the X = 0 plane,' Y = 0 plane, and 

Z - 0 plane, and the model simulated a vessel with two symmetric surface cracks (180 degrees 
apart). The X = b plane was free. 

The stress-intensity factor solutions were obtained by solving the complementary problem 
of applied stresses on the crack surfaces. Four applied stress distributions on the crack 
surfaces were analyzed: constant, linear, quadratic, and cubic. These stresses, which were 

applied to the crack surfaces as shown in figure 2(c), were symmetric about the y = 0 plane 
and were given by 


’2 ■ (if 


for j = 0 to 3 


( 1 ) 


where z is measured from the crack mouth toward the crack front. Solutions for these four 
stress distributions were superimposed to obtain stress-intensity factors for the pressurized 
cylinder. (These four solutions can be superimposed to obtain stress-intensity factors for 
other stress distributions, such as those caused by thermal shock.) 
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4. Stress-Intensity Factor 


The Mode I stress-intensity factor, Kj, at any point along the surface crack was taken 

to be 



for j - 0 to 3. Gj is the boundary-correction factor corresponding to the J th stress 
distribution from eq. (1). Q, the shape factor for an elliptical crack, is given by the 
square of the complete elliptic integral of the second kind. The vessel length (2b) was 

always chosen large enough that the length would have a negligible effect on stress intensity 
(b/c >10). 

The stress-intensity factors for a surface crack in an internally pressurized cylinder 
were obtained by appropriate superposition of the results given by eq. (2). For convenience, 
the stress-intensity factor was written as 


£R 

t 


'TTr F 




(3) 


where pR/t is the hoop stress and F is the boundary-correction factor for a surface crack 
on the inside of an internally pressurized cylinder. The expression for F, in terms of G , 
was obtained from the first four terms of a power-series expansion of Lame's solution [10] 


for the hoop stress in an internally pressurized cylinder. The result is 



where each was obtained from the appropriate finite-element solution. The correction 

factor F includes also the influence of the internal pressure, p, acting on the crack sur- 
faces. Values for F were calculated as a function of a/c, a/t, and <f> for t/R values 
of 0.1 and 0.25* The a/c ratios were 0.2. 0.4, and 1; and the a/t ratios were 0.2, 0.5, 
and 0.8. 

The stress-intensity factors from the finite-element models were obtained by using a 
nodal-force method, details of which are £iven in references [7] and [9] . In this method, 
the nodal forces normal to the crack plane and ahead of the crack front are used to evaluate 
the stress-intensity factors. 


5. Results and Discussion 

In the following sections, results are presented for two symmetric semi-elliptical sur- 
face cracks on the inside of pressurized cylinders. The stress-intensity factor variations 
along the crack front for various surface cracks (a/c =0.2 and 1) are presented as a func- 
tion of a/t. An empirical equation for the stress-intensity factor is also developed for a 
wide range of configuration parameters. An estimate for the stress-intensity factor for a 
single surface crack is also presented. The stress-intensity factors are compared with other 
solutions from the literature. 

5. 1 Semi-Circular Surface Cracks 

Figure 3 shows the boundary-correction factors for two symmetric sesni-circular surface 
cracks (a/c * 1) as a function of the parametric angle, <J>, and a/t for t/R = 0, 0.1, and 
0.25. For t/R * 0 (flat plate [8,9]), the pR/t stress in eq. (3) is replaced by S t , a 
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remote uniform applied stress. For a fixed t/R ratio, the correction factors are higher 
for larger a/t ratios. Also, for a fixed a/t ratio, the correction factors are higher 
for smaller t/R ratios. The maximum correction factor (or stress-intensity factor) 
occurred at the intersection of the crack front with the inner surface (<{> *= 0) . 

5.2 Semi-Elliptical Surface Cracks 

Figure 4 shows the boundary-correction factors for two symmetric semi-elliptical surface 
cracks (a/c = 0.2) as a function of * and a/t, for t/R - 0, 0.1, and 0.25. Again, for a 
fixed t/R ratio, the correction factors are higher for larger a/t ratios. For a <given 
a/t ratio, smaller t/R ratios gave higher correction factors. In contrast to results for 
a semi-circular surface crack, the maximum correction factor (or stress-intensity factor) 
occurred at the maximum depth point (<£ - tt/2) . 

5.3 Stress-Intensity Factor Equation 

The results shown in figures 3 and 4 suggest that the ratio of the correction factors 
for a given t/R and those for a flat plate (t/R = 0) are nearly independent of the para- 
metric angle, so that the curve for t/R = 0 can be scaled to approximate the results for 
t/R =0.1 and 0.25. Figure 5 shows f £ (the ratio of F for a given t/R to F Q ) as a 
function of t/R. F Q is the correction factor for a flat plate [8,9]. The nbols (and 
bars) give the average (and range) of f for a given valee of a/t with a/c = 0.2, 0.4, 
or 1 and any value of <f». (For clarity, results for a/t = 0.5 are not shown.) These 
results were found to be closely approximated by 


f 

c 



+ R 



•f 


1 



(5) 


In figure 5, the upper curve shows the exact limiting solution for a/c = 0 and a/t *= 0. 
The upper curve was obtained from Lame's stresses [10] on the inside of the cylinder and the 
solution for an edge-crack in a semi-infinite plate, and is given by eq. (5) with a/t - 0. 
The other curves show results from eq. (5) for various a/t ratios. 

The stress-intensity factor for two symmetric surface cracks on the inside of a pres- 
surized cylinder is given by eq. (3) where the following approximate expression for F has 
been fitted to the present results and those of references [8] and [9] : 


• 9, [ m i + ^(f ) 2 ♦ M 3 (?) 3 ]f' 


F = 0 


M. = 1.13 - 0.09 - 
1 c 


pf A f 
$ c 


( 6 ) 

(7) 


M 2 = 0.54 + 


0.89 

0.2 + ® 

c 


M 3 = 0.5 - 


0.65 + ~ 
c 


+ 14 


(■ - 7) 


24 


( 8 ) 


(9) 


1 + [0, + 0.35(|) 


(I - sin <$Y 


g 


( 10 ) 


and f . 
v 



( 11 ) 


for 0<^<0.8, 0 < ^ < 1, t/R < 0.25, and any <{>. Equation (6) is within about ±5 per- 

cent of the finite-element results. 

An estimate for the stress- intensity factors for a single crack in terms of two sym- 
metric cracks was obtained from an equation (eq. (2)) given in reference [11]. For a/c * 0 
and t/R * 0.1, the stress- intensity factor for a single crack was about 2 percent lower than 
that for two symmetric cracks, whereas for t/R = 0.25 the stress- intensity factor for a 
single crack was about 4 percent lower than that for two cracks. For larger a/c ratios the 
differences are smaller. Thus, eq. (6) can also be used for a single surface crack without 
appreciable error. 

5. 4 Comparisons with Other Solutions 

Figure 6 shows the stress- intensity boundary-correction factors obtained by several 
investigators for a semi-elliptical surface crack (a/c = 1/3 and a/t s 0.8) on the inside 
of a cylindrical vessel with t/R =0.1. The crack surfaces were subjected to stress dis- 
tributions given by eq. (1). The correction factors, H. , are those used in references [5] 
and [6] , and are related to by 

B G 

H *= - ^ for j = 0 to 3 02) 

where B Q = 1, B J = tt/2, B 2 * 2, B 3 * 3tt/ 4, and f is given by eq. (11). The present 
results are shown as symbols. The results from reference [6], obtained from the boundary- 
integral equation method, are shown as solid curves and are in good agreement (±2 percent) 
with the present results for <{> > tt/4. The dashed curves show the results from reference [5]. 
Their results were obtained from the finite-element method and were within ±S percent of the 
present results. (For clarity, the results for were not shown.) 

6. Concluding Remarks 

A three-dimensional finite-element elastic stress analysis was used to calculate stress- 
intensity factors for a wide range of semi-elliptical surface cracks on the inside of 
cylindrical pressure vessels. The ratio of crack depth to crack length ranged from 0.2 to 1; 
the ratio of crack depth to wall thickness ranged from 0.2 to 0.8; and the ratio of wall 
thickness to vessel radius was 0.1 and 0.25. Singularity elements were used along the crack 
front and linear-strain elements were used elsewhere. The models of these configurations had 
about 6500 degrees of freedom. A nodal-force method was used to evaluate the stress-intensity 
factors. 

The stress-intensity factors for surface cracks in pressurized cylinders were similar to 
those calculated for surface cracks in flat plates as a function of the parametric angle. 

For scmi-circular cracks, the stress-intensity factors were maximum at the Intersection of the 
crack with the inside surface of the cylinder, but for semi-elliptical cracks the values were 
largest at the maximum depth point. Larger crack-depth-to-wall-thickness ratios and larger 
wall-thickness-to-vesse 1-radius ratios gave higher stress-intensity factors for all surface- 
crack configurations considered. 
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The present results were compared to other analyses of internal surface cracks in 
cylinders. The surface-crack configuration had a crack-depth-to-vall-thickness ratio of 0.8, 
a crack-dcpth-to-crack-length ratio of 1/3, and a vall-thickness-to-vessel-radius ratio of 
0.1. The cracks were subjected to constant, linear, quadratic, and cubic stress distribution 
on the crack surfaces. The results from a boundary-integral equation method were in good 
agreement (generally ±2 percent) and those from a finite-element method were in fair agree— 
ment (±8 percent) with the present results. 

An empirical equation for the stress-intensity factors for an internal surface crack in 
a pressurized cylinder was developed to estimate the present results. The equation applies 
over a wide range of conf iguration parameters and was within about ±5 percent of the present 
results. The results obtained herein should be useful in correlating fatigue-crack-grovth 
rates and in calculating fracture toughness for the surface crack in a pressurized cylinder. 
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